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\ We examine the decay modes rj/rf — ► 7r + 7r~7 within the context of the Hidden 

Local Symmetry (HLS) Model. Using numerical information derived in previous 
fits to VPj and Ve + e~ decay modes in isolation and the p lineshape determined 
in a previous fit to the pion form factor, we show that all aspects of these decays 
can be predicted with fair accuracy. Freeing some parameters does not improve 
the picture. This is interpreted as a strong evidence in favor of the box anomaly 
in the 77/7/ decays, which occurs at precisely the level expected. We also construct 
the set of equations defining the amplitudes for rj/rj' — * tt + tt~j and 77/77' — ► 77 at 
^ ■ the chiral limit, as predicted from the anomalous HLS Lagrangian appropriately 

broken. This provides a set of four equations depending on only one parameter, 
instead of three for the traditional set. This is also shown to match the (two-angle, 
two-decay-constant) 77 — 7/ mixing scheme recently proposed and is also fairly well 
fulfilled by the data. The information returned from fits also matches expectations 
from previously published fits to the UP7 decay modes in isolation. 
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1 Introduction 



Interactions and decays of light mesons fit well within the framework of Chiral Pertur- 
bation Theory (ChPT) pQ. Strictly speaking, the ChPT framework applies to the octet 
members of the pseudoscalar sector (n,K,r)s) which behave as Goldstone bosons whose 
masses vanish at the chiral limit. Relying on the large N c limit of QCD, an extended 
ChPT framework (EChPT) has been defined [2HH| including the singlet 770 state which 
keeps a non-zero mass at the chiral limit, but this vanishes in the large N c limit. On the 
other hand, the decays 7r°/rj/rj' — > 77, are understood as proceeding from the so-called 
triangle anomaly. These are accounted for by means of the Wess-Zumino-Witten (WZW) 
Lagrangian [3J|S] which is also normally incorporated into the ChPT Lagrangian j2HH]- 

Other anomalous processes describing the (7r°7r + 7r~7) vertex and the decay mode (rj — > 
7r + 7r~7) have been identified long ago within the context of Current Algebra t 6j ; they are 
presently referred to as box anomalies. Triangle and box anomalies are now derived from 
the WZW Lagrangian. The box anomaly part of the WZW Lagrangian predicts exactly 
the values of the amplitudes for the couplings (7r°7r + 7r~7), (7]n + n^j) and (7/7r + 7r~7) at 
the chiral limit ; however, the momentum dependence of the corresponding amplitudes 
is not predicted and should be modelled. When dealing with experimental data, this 
momentum dependence is naturally accounted for by vector meson contributions and, 
then, the question becomes whether these alone account for the box anomalies or whether 
an additional contact term (possibly simulating high mass resonances) is needed ; if this 
contact term (CT) is needed, it should have a definite value in order to stay consistent 
with the rigorous predictions of the WZW Lagrangian. 

Therefore, from an experimental point of view, the question of the relevance of the 
box anomaly phenomenon turns out to check the need for a well-defined contact term 
besides the usual resonant contributions. This question is still awaiting a definite and 
unambiguous signature. 

In its simplest figure, the problem of the relevance of the box anomaly phenomenon 
is adressed in the coupling (7r°7r + 7r~7). The relevance of a possible contact term beside 
vector meson exchanges has been examined. A value for this coupling has been extracted 
from experimental data [7] and found close to expectations 1 (only 2a apart). 

A cleaner environment could be provided by the decay modes Tj/rf —> 7r + 7r~7 which 
are also accounted for in the WZW Lagrangian. Several pieces of information are avail- 
able : the partial widths [§| are known with an accuracy of the order 10% cross-checked 
by several means, the 77 spectrum as function of the photon momentum has been mea- 
sured long ago fOm] an d provides useful information. Finally, measurements of the 
rf spectrum as function of the dipion invariant-mass have been performed twelve times, 
with various levels of precision, and the corresponding data are published as papers from 
several Collaborations [13111111311311111111 or are available as PhD theses H3ED1EI1E2I • 
The latest measurements have been performed recently by CERN Collaborations [231l24| . 

This already represents a large amount of information covering all aspects of these 
decays. This should allow a reasonably well founded analysis in a search for the box 

1 See, also, the discussion in 
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anomaly in the rj/rf system. 

As stated above, predictions on box anomalies are given at the chiral point and 77/77' 
spectra clearly extend in regions where accounting for the p exchange cannot be avoided 
in order to match experimental information. The magnitude of the f]/r]' — > 7r + 7r~7 
partial widths should also be influenced by the p exchange. Stated otherwise, including 
momentum (invariant-mass) dependence within their spectra is essential and should be 
done in a consistent way in order to extract reliably information on the box anomaly 
from data. 

Indeed, in these decays, two contributions are a priori competing : a contact term 
and a (dominant) resonant one -the p exchange- with contributions of (sometimes) very 
different magnitudes. In order to detect reliably the former, the latter has to be known 
with enough accuracy and, possibly, fixed. The sharing in the anomalous amplitudes at 
the chiral limit between the contact term and the resonant term might also have to be 
understood unambiguously. 

Therefore, a global framework is needed where the vector meson degrees of freedom 
are explicitly accounted for together with pseudoscalar mesons and the contact terms. 
Several such frameworks implementing vector dominance (VMD) in effective Lagrangians 
have been defined : Resonance Chiral Perturbation Theory |HI2H], massive Yang-Mills 
fields [2111211121], Hidden Local Symmetry (HLS) Model [2HH20]- It was soon shown [21] 
that all these approaches were physically equivalent. For convenience, we work within 
the HLS Model context. 

A second issue makes the difference between the tt° box anomaly and the 77/77' ones. 
The former is practically insensitive to symmetry breaking effects (Isospin Symmetry 
breaking is a small effect), the latter however sharply depends on how SU(3) Symmetry 
and Nonet Symmetry breakings really take place. Therefore a reliable breaking scheme 
of the 77/7/ sector should also be defined and checked in the triangle and box anomaly 
sectors. It should also be validated in all processes where it has to apply, like V — > P7 
and P —>■ V'-f decays. One has already noted some confusion [22] in the meaning of the 
decay constants entering the amplitudes for the rj and rf decays to two photons. 

If one limits oneself to collecting some VMD term for the p contribution (even if 
motivated) and simply adds it with a phase space term to be fit, one can be led to 
ambiguities [221121] when solving the Chanowitz equations [22] which represent the tra- 
ditional way of describing the 77/77' mixing (see also [36J). Along the same line, if the 
breaking scheme generally used j3lIll3l)ll3l)ll3T|l3*5] happens to be inappropriate in order to 
describe the 77/77' system, extracting the box anomaly constant values from data becomes 
hazardous. 

A scheme for implementing SU(3) symmetry breaking in the full HLS Lagrangian 
has been already defined [23110] ■ This scheme, referred to as BKY has been proved [22] 
to meet all (E)ChPT requirements and allows a successful account of a very large set of 
experimental data [HP^] ■ A brief global account of the full breaking scheme we advocate 
is summarized in Appendix A to [12] ■ The non-anomalous sector has been used in pion 
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form factor studies providing also consistent results [HJII3]. 

Therefore, in this paper, we intend to extend the realm of the broken HLS model by 
studying the decays 77/77' — > 7r + 7r~7. The behaviour of the model can then be examined 
in a context where the box anomaly phenomenon is expected to be present. One can 
hope extracting unambiguously the information about the relevance of this phenomenon 
from experimental data. 

The paper is organized as follows ; in Section |21 we recall the traditional expressions 
of the decay amplitudes at the chiral limit for rj/rj' — * 77 and 77/7/ — > 7r + 7r~7. In 
Section EJ we outline the derivation of the full anomalous sector of the HLS Model, 
mostly refering to the basic paper [30! • ^ n Section 0] we recall shortly how the BKY 
breaking of SU(3) Symmetry and the breaking of Nonet Symmetry has been performed 
and tested. In Section we recall the result of applying this to 77/77' — > 77 as it provides 
an unconventional set of expressions for the amplitudes at the chiral limit. 

In Section we develop the structure predicted for the rj/rj' — > 7r + 7r~7 decay modes 
by the broken HLS Model. We first show that the BKY breaking scheme provides also 
unconventional expressions for the box anomaly amplitudes at the chiral limit. We also 
show that all information related with these decay modes (parameters and p meson 
lineshape) have been already derived numerically and functionally in other sectors of the 
low energy phenomenology. It thus follows that all properties of the 77/77' — > 7r + 7r~7 decay 
modes can be predicted without any numerical or functional freedom. In Section we 
examine the predictions of this model for the 77/77' — > 7r + 7r _ 7 partial widths and for their 
dipion invariant-mass spectra. 

After reviewing shortly the status of the available experimental data on this subject in 
Section |Hl we devote Section |H1 to comparing the predicted lineshapes with the published 
experimental spectra. Section [TU] is devoted to performing a global fit of the shape and 
yield information for the rj/rj' — > 7r + 7r~7 modes in order to check precisely the relevance 
of the numerical parameters which were all fixed from analysis of other independent data 
sets. In Section [TTJ we propose, for comparison, fits of the anomalous amplitudes at the 
chiral limit, under various conditions and show that the one (instead of three, usually) 
parameter dependence of these gets a strong support from data. 

Finally, Section^Jis devoted to a summary of the results obtained and to conclusions. 

2 Radiative Decays of Neutral Pseudoscalar Mesons 

Some interactions (or decay modes) of neutral pseudoscalar mesons (P = 7r°,77, 77') 
are described by matrix elements having the wrong parity and are called anomalous. 
Anomalous interactions were treated by Wess and Zumino j3] and then expounded upon 
by Witten [3] ; they are given by the anomalous action, which we shall refer to as Twzw- 
For the purpose of this paper, two pieces 2 from r WZ w are relevant : 

2 Here and in the following, we denote by V the (massive) vector field matrix, by A the electromagnetic 
field and by P the pseudoscalar field matrix. The matrix normalization we use for these have defined 
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with e 2 = 47ra, and f n = 92.42 MeV ; Q is the quark charge matrix given by 3 A is the 
electromagnetic field and P is the bare pseudoscalar field matrix. From there, amplitude 
intensities at the chiral point can be derived. 

The first piece £ 77 p describes the decays 7T°/r]/r]' — > 77. The second piece £ 7 ppp 
contains an interaction term 7 — > 7r + 7r°7r~ briefly discussed in the Introduction. This last 
piece contains also terms which account for the anomalous decay modes 77/77' — > 7r + 7r~7. 

Without introducing symmetry breaking effects, the Lagrangian pieces in Eq. (0) 
can give reliable predictions for processes involving only pions. In order to deal with 
interactions involving 77 or 7/ mesons, one has to feed these Lagrangians with SU(3) 
and Nonet Symmetry breaking mechanisms. Usually these breaking mechanisms are 
considered to arise from the naive replacement of the pseudoscalar decay constants [SSJ 
OHIIHIEZIIIH] • Using obvious notations, the amplitudes at the chiral point derived from 
Eqs. ((*[)) can be written : 



T{X - 77) = B x (0) e^e^kpk'p 

T (X -> 7T+7T-7) = E x (0) e^e p k u p+p- 
(X = 77,7/), where the coefficients are, assuming N c = 3 : 
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using the traditional one-angle mixing scheme. The procedure is thus obvious : one 
replaces one power of by the octet (/ 8 ) or singlet (/ ) decay constant understood 
under their customary definitions in (Extended) ChPT. In the following, we refer to 
X — > 77 and X — > 7r + 7r~7 as triangle and box anomalies. 

in |29ll4()ir^] ; our normalization for the SU(3) flavor matrices differs from those in [H] by a factor of 
2 : ^Hoistein = Twitten = 2 Tfj LS . Moreover, we use without distinction VPj and AVP to name the 
corresponding coupling. 

3 There is an intimate connection between the charge of quarks and the value of N c in the anomalous 
action gSlEHET] ; Q = Diag(2/3 -1/3,-1/3) if N c = 3.. 
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This implies several assumptions which are traditionally made in an implicit way 
0133 : 



• the decay constant /§ and /o are the (usual) decay constants of ChPT defined from 
current expectation values : (0\J^\r] s (q)) = ifsq^, and (0| J®\rio(q)) = if^q^, 

• SU(3) and Nonet Symmetry Breakings act in exactly the same way for the triangle 
and box anomalies. 

These equations have been used in several ways and they underlie decades of phe- 
nomenological work on the 77/7/ mixing. For instance, Refs. |3EIIHIISZ| consider the two- 
photon decay widths of the 77 and 7/ mesons and the ratio fs/ f n — 1-3 from ChPT £Q to 
derive 9p ~ —20° and f / f w ~ 1.04 ; this meets ChPT expectations [2HH| if one identifies 
9 P with the presently named 9 8 . Comparable results j^HIHE] are derived by using the 
four equations above, after extracting the box anomaly constants from the dipion mass 
spectra in 77/7/ — > 7r + 7r~7 decays. The third Eq. (jHJ) has also been used with accepted 
parameter values (/ 8 / fn = 1-25, f / f w = 1.04 and dp = —20.6°) inside the HLS Model 
to derive a successfull description of 77 — ► 7r + 7r~7 in isolation [15] . 

The validity of the first two Eqs. (J3J) has been recently addressed and consistency of 
these with the 77/77' breaking scheme derived from EChPT [21IH] has been found doubtfull 

There is no currently known examination of the last two Eqs. (J3J) ; however, some 
remarks on the renormalization of the WZW box term |40] tend to indicate that these 
are also doubtful. Therefore, the phenomenological results derived from using Eqs. © 
have to be reexamined in a consistent framework. 

3 The Anomalous Sector in the HLS Model 

The HLS Model originally describes the 7 — V transitions, all couplings of the kind 
VPP and possibly APP, if the specific parameter a of the HLS model [2D., is not fixed 
in order to recover the traditional VMD formulation (a = 2). In this framework, the 
main decay mode u — > 7r + 7r°7r~ of the uj meson is, for instance, absent as clear from the 
explicit expression of the HLS Lagrangian |4T)] . 

As seen above, anomalous interactions involving pseudoscalar mesons and photons 
are contained in Twzw [HE]- These terms were included in the Hidden Local Symmetry 
Lagrangian by Fujiwara et al. [3U], along with anomalous vector meson (V) interactions in 
such a way that the low energy anomalous processes (in the chiral limit where m n = 0) 
7 — > 7r + 7r°7r~ and 7r° — > 77 are solely predicted by Twzw- The construction of this 
HLS anomalous Lagrangian, originally performed in [HOI ' ^ s discussed in detail in several 
excellent reviews jSHUSH]- Here, we limit ourselves to a brief outline of its derivation, 
pointing out the motivation for some important assumptions. The anomalous action has 
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the form : 

r = r wzw + Tfktuy 

(4) 

Tfktuy = J2i=i c % J d x Ci 

where the Cj are entirely arbitrary constants. The Lagrangians Ci are given in [30j and 
each of them contains APPP and AAP pieces which would contribute to the anomalous 
decays, but are cancelled by APV terms. These Lagrangians contain also VPPP and 
VVP pieces [303 - 

In view of extending the assumption of dominance of vector mesons (VMD) to the 
anomalous sector, it was first shown that a set of Cj in Eq. can be defined in such 
a way that tc° — > 77 occurs solely through tt° — > ujp° followed by the (non-anomalous) 
transitions uj — > 7 and p° — ► 7. The tt° width thus derived is identical to the Current 
Algebra prediction reproduced by £ 77 p defined in the previous Section. 

It was also shown [HOj that complete vector dominance can be achieved, where the 
direct term APPP is converted to VPPP, with some other set of q. In this framework, 
the decay u — > 7r + 7r°7r~ occurs through uj — > 7r°p° (the VVP term) followed by p° — > 7r + 7r~ 
and through the contact term (VPPP) which gives a direct contribution uj — > 7r + 7r°7r~. 
However, it happens that the VVP contribution (which provides alone the correct u 
partial width) is numerically reduced in a significant way by the contact term (VPPP). 
In view of this, j3U] proposes another set of q which provides an anomalous effective 
Lagrangian containing only a VVP term and, besides, the standard WZW term APPP 
in the following combination : 

£ FKUTY = -^^TtlWAVrFl ~ . (5) 

47T Z Z 

where £ 7 ppp is defined in Eqs. (0). One should note that the normalization affecting the 
WZW part of this Lagrangian is a pure prediction of the HLS Model based on a definite 
extension of the VMD concept to anomalous processes. 

Focussing on decays like 77/7/ — > 7r + 7r - 7, one readily sees from this expression that, 
in order to recover the behaviour expected from £ 7 ppp alone, these two terms should 
contribute to the box anomaly (i.e. the full amplitude at the chiral limit) in the following 
ratio : 

VMD : CT = -3:1 , 

at the chiral limit ; "VMD" names here the contribution generated by the first term 
in Eq. (jSJ) and "CT" (contact term) those generated from the second term. Thus, the 
"VMD" contribution, generated by the triangle anomaly generalized to VVP couplings, 
is predicted dominant at the chiral limit. 

Within this framework, the main uj decay mode proceeds only from uj —>■ p°7r° followed 
by p° — ► 7r + 7r~ and <p — ► tt + 7t~7t° proceeds solely from uj — <ft mixing. The experimental 
situation concerning the decay mode — > 7r + 7r~7r° is conflicting. Indeed, a recent result 
from the SND Collaboration claims for no significant evidence in favor of a contact 
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(f) — > 7T + 7r _ 7r term in their e + e~ — > 7r + 7r~7r° data and provides a new upper bound 
much more stringent than previous ones [H] ; however, using their own data on the same 
physical process, the KLOE Collaboration [SS] claims that a significant contact term is 
present in their data. Actually, as there is currently no available analysis performed using 
consistenly a full VVP and VPPP Lagrangian or a Lagrangian like in Eq. (JHJ), no founded 
conclusion can really be drawn. 

Processes like 7r°/r)/r)' — > 77 occur solely through n /r]/r]' — ► W followed by V, V — > 
7. However, transitions like 7 — > 7r + 7r°7r~ or decays like 77/77' — ► 7r + 7r~7 receive contri- 
butions from the contact term and from the VVP term (essentially through p meson 
exchange) . 

The VVP piece of this effective Lagrangian has been used successfully in several recent 
studies (H1I12IE2IES1 and proved to predict (after implementing appropriate symmetry 
breaking mechanisms) up to 26 physics information with a number of free independent 
parameters ranging from 6 to 9 (when Isospin Symmetry breaking is considered [53]). 

4 The Extended BKY Symmetry Breaking Scheme 

The study of SU(3) breaking of the HLS Model has been initiated by BKY [39J who 
proposed a simple and elegant mechanism. However, the procedure was soon understood 
as breaking also Hermiticity of the derived Lagrangian, which was clearly an undesired 
feature. A slight modification [3U] of the original BKY procedure was shown to cure this 
problem and to produce a quite acceptable broken Lagrangian (see Eq. (A5) in |40j). 
The way field renormalization has to be performed turns out to define the renormalized 
field matrix (denoted P') in terms of the bare field matrix (denoted P) by : 

P = X A 1/2 P'X A 1/2 , (6) 

where the breaking matrix is Xa = Diag(l, 1, z), with z = [fx/ fir} 2 - 

As such, the (original) BKY breaking scheme can only address a limited amount of 
physics processes, as all information related with the 77 meson can only be treated crudely 
and the properties of the 77' meson cannot be addressed. 

In order to address physics information about the 77/77' system appropriately, the 
singlet sector has first to be introduced in the original HLS Lagrangian. This has been 
done by using jlHI the U(3) symmetric field matrix P = Pg + Pq instead of only P = P$ 
when constructing the Lagrangian. This is found to provide the HLS Lagrangian with 
the canonical kinetic energy of the singlet state (770 field) while this does not modify the 
interaction Lagrangian [3U] by adding 770-dependent pieces. 

The step further is to break the Ua(1) symmetry by introducing determinant terms 
[51)] into the effective Lagrangian which becomes jH2] : 

£ = £hls - \^ 2 vl + ttA^t/q^t/o (7) 
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By means of this (modified) BKY breaking scheme, the HLS Lagrangian can now ad- 
dress the r]/rf system with a complete scheme of symmetry breaking (SU(3) and Nonet 
Symmetries). A 770 mass term is generated and the kinetic energy term of the Lagrangian 
is modified in a non-canonical way, which implies that a field transformation to renormal- 
ized fields has to be performed. This can be done through the two-step renormalization 
procedure defined in [32] and outlined in the Appendix. This transformation is well 
approximated [32] by : 

P = X A ll2 lpi + xP' }X A 112 ■ (8) 

this has been shown to differ [32] from the exact field transformation only by terms of 
second order in the breaking parameters (\z — 1], [x — 1]). 

This transformation 4 was postulated (or fortunately anticipated) in 0T| in order 
to study the full set of AVP radiative decays, especially the modes involving the rj/rj' 
mesons. Using this transformation has provided a fairly good description of the whole 
physics accessible to this broken Lagrangian [1IJ[32J[32JES] with only a small number of 
parameters, as already noted. 

Departures from Nonet Symmetry were observed |H] by extracting from data x = 
0.90 ± 0.02, significantly different from 1. One may rise the question of the correspon- 
dance between x in Eq. (JHJ) and the basic Nonet Symmetry breaking parameter A of the 
Lagrangian in Eq. (JJJ) - which is [12] nothing but the Ai parameter of [2HH]- One can 
write : 

Indeed, as shown in the Appendix, x absorbs a small influence of SU(3) symmetry 
breaking (about 5% of its fitted value). 



Comparison of all results of this broken HLS Lagrangian, especially decay constants 
and mixing angles, with the available (E)ChPT estimates of the same parameters |2*|I3|I5U] 
was done and appeared also fully satisfactory. It is worth remarking that the HLS 
phenomenology was yielding an estimate for the (E)ChPT mixing angle #0 much smaller 
in magnitude than the (E)ChPT leading order estimate (—0.05° ± 0.99° in contrast with 
~ —4°), but in fair correspondence with a more recent EChPT next-to-leading order 
calculation [58J which yields 6 = [-2.5°, +0.5°]. 

It is also worth remarking that the (full) breaking scheme just outined anticipated [1J 
the branching fraction for — > 7/7 with a value twice smaller than its contemporary 
measurement [53] . This predicted value coincides with all recent measurements performed 
with larger statistics 0. 



The quasi-vanishing of #0 has two interesting consequences. On the one hand, it 
allows to relate the traditional wave-function mixing angle with the recently defined #8 

4 The motivation behind this postulate was that weighting differently the singlet and octet parts of 
the P and V field matrices allows to derive the most general parametrization J37] of the VPj cou- 
pling constants consistent with only SU(3) symmetry in the vector and pseudoscalar sectors, while the 
corresponding U(3) symmetries are both broken. 
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mixing angle by providing ~ 29 p (fulfilled at a few percent level) ; the derivation 
is given in the Appendix for the exact field transformation. 

On the other hand, the condition Qq — relates the Nonet symmetry breaking pa- 
rameter x to Op : 

(1 — z) 

tanflp = V2 2 + J x (10) 

This relation is fulfilled with a high degree of numerical accuracy [32] and only reflects 
that the ChPT mixing angle is not significantly affected by symmetry breaking effects. 
This relation will be somewhat refined (See Section ^2 and the Appendix). 

It then follows that, from the three originally free breaking parameters associated 
with the pseudoscalar sector (z, x, Op), only one remains unconstrained. It could be 
either of x or Op ; however, it will be shown that x might be prefered. 

We do not go on discussing here symmetry breaking in the vector meson sector as it 
is not in the stream of the present paper ; we refer interested readers to 1^1 ) 112] 
where this is discussed in details. 

Some remarks are of relevance. The combined Nonet Symmetry and SU(3) breaking 
scheme of the HLS Lagrangian presented in this Section defines what we name the Ex- 
tended BKY breaking scheme. It restores the relevance of a one angle mixing scheme 
for the rj/rf system. However, this does not give any support to the traditional break- 
ing scheme as expressed by Eqs. (J3J). In contrast, it matches fairly well all expec- 
tations of the two-angle, two-coupling constant mixing scheme recently derived from 
EChPT OEJEDHin] at leading order in the breaking parameters. 

This full breaking scheme is also mathematically equivalent to the recently proposed 
|5U1Iq"T] breaking in the quark flavor basis framework ; it might be prefered as a definite 
concept like Nonet Symmetry breaking, which underlies some Lagrangian pieces (£2) 
of EChPT, can be implemented more clearly and tracing its effect in phenomenology is 
easier. 



5 Two— Photon Decay Widths of the rj and rj mesons 

The two-photon decay widths of the rj and 7/ mesons can be derived easily from 
the HLS Lagrangian (the VVP part of Eq. (jSJ)) supplemented by the V7 transition 
amplitudes of the non-anomalous HLS Lagrangian) after renormalizing to physical fields 
by Eq. (JSJ). Applying directly the same Eq. (jSJ) to the WZW Lagrangian £ 77 p in Eq. 
(JTJ) leads exactly to the same result 5 j32] : 



0,(0) 
GV(0) 



a 



a: 



5z 



71 



3z 
5z-2 
3z 



5z + 1 

cos Op - V2 — - — x sin Op 
2>z 



5z + 1 

sin — x cos Op 

3z 



IV. 



5 From now on, we assume N c = 3. 
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These expressions compare well with the corresponding quantities in Eqs. (jSJ)- How- 
ever, this correspondence is only formal as, defining / 8 and f by : 

U 5z-2 f n 5z + l 

= — j = — x > (12) 

h 3z /o Qz 

yields / 8 = 0.82 f n (and f = 1.17 f n ), which has little to do with numerical expectations 
from ChPT (extended or not). It was proved in that these are not the standard 
EChPT decay constants. These can be derived from our broken Lagrangian, yielding 
information which matches jH2] fairly well EChPT expectations 

This proved the basic consistency of breaking scheme presented in the previous Section 
with EChPT. The formulation given in Eq. (fTT|) could look, at leading order, more 
tractable than present standard expressions ; it makes indeed much clearer that the 
number of parameters to be determined phenomenologically is limited. 

Therefore, the first basic assumption which underlies the understanding of Eqs. (jSJ) 
is not fulfilled by the BKY breaking scheme p^HHTi] and this is independent of whether 
Nonet Symmetry is broken. 

Let us remind that Eqs. (fTTjl give the two-photon radiative decay widths of the 77/7/ 
mesons with good accuracy. These can even be predicted by using solely the value of x 
extracted from fit to the (independent) set of AVP decay modes of light mesons |32j . 
Fixing z = [fx/ fir} 2 to its experimental value and assuming Eq. (JTUJ), Eqs. (fTTj) become a 
constrained system depending on only one parameter and can be solved providing results 
consistent with using, instead, the AVP decay mode information. 

Of course, the mixing angle dp entering Eqs. does not coincide with 8g and 

is derived as 9p = —10.32° ± 0.20° when requiring the constraint Eq. (fTU|) to hold 
exactly ; the corresponding value for 6$ ~ —20° compares well to expectations fTll2l'jj . One 
should note a recent estimate of 9 P = —10° ± 2° provided by lattice QCD computations 
which strongly supports this phenomenologically extracted value. 

Therefore, the picture represented by Eqs. (|llj). which does not meet traditional 
expectations [H1EHIEHIE3] > matches quite well all relevant information from ChPT and 
QCD, and, last but not least, corresponds to a satisfactory description of the whole set 
of two-body radiative decays of light mesons [rljU^l EJ^] . 



6 The HLS Model For 77/77' tt+tt ^ Decay Modes 

Using the effective Lagrangian in Eq. (jSJ), the processes 77/77' — > 7r + 7r _ 7 receive VMD 
contributions from the VVP term and CT contributions from the C^ppp piece. The 
purpose of this Section is to examine carefully these decay modes. These will also lead 
us to question the last two (box) anomaly equations Eqs. ©• 
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6.1 Basic Lagrangians 



Within the HLS Model, the VVP part of rj/rj' — > 7r + 7r~7 involves, beside the p meson, 
the interplay of the u and <fi mesons to their decay mode to titi only 6 . However, these 
(isospin violating) couplings are small enough to be safely neglected. Additionally, the (f) 
meson is outside the decay phase space of both 77 and rf mesons, and the accuracy of the 
data is far from allowing any uj effect to be significant or simply visible in the 7/ dipion 
invariant-mass spectrum. 

It can be shown [12], that the VP'y couplings following from the anomalous sector of 
HLS model can be derived from the corresponding (VP7) piece of : 



£ = Ce^TT[d,(eQA u + gV u )d p {eQA a + gV a )X A 1/2 (P^ + xP^)X A 



-1/2, 



(13) 



where g is the universal vector coupling of the HLS Model [29]. The value for C = 
— 3/(47r 2 / 7r ) is fixed by normalizing the AAP term in Eq. (|13[) to the corresponding 
WZW Lagrangian in Eq. |T|). 

This equation could essentially be considered as a way to postulate VMD for VVP 
couplings and it also gives the normalization shown in Eq. (jlj]). Focussing on the piece 
of Eq. (fT3*|) related with neutral pseudoscalar mesons, one gets : 



C 



7pP° 



eg 



1 ^3 13 
2 n + ^ r]8 + l2 m 



(14) 



with obvious notations. The p meson decay relevant for the present study is driven by : 



r _ 



7T <9 M 7T + — TT+d^TT' 



(15) 



which can be extracted from the non-anomalous (broken or not) HLS Lagrangian. The 
parameter a is the specific HLS parameter expected such that a = 2 in traditional 
formulations of VMD, but has always been fitted in the range a = 2.3 4- 2.5 from radiative 
and leptonic decays of light mesons [5TJI32] and in pion form factor studies ^SHUO] ■ ^ e 
recall that the p mass in the HLS Model is not free and is given by the (extended) KSFR 
relation m 2 p = ag 2 f 2 which does not coincide with traditional mass values (~ 830 MeV 
versus ~ 775 MeV) ; these happen to be only a matter of definition [32] ■ Finally, the 
CT contributions are contained in the following Lagrangian piece extracted from Eq. (0) 
-for the normalization- and Eq. (JTJ : 



C 



- P o 



[——e^A^d^ 



(16) 



6 Indeed, the non-anomalous HLS Lagrangian, broken or not ^U], contains no couplings like rjirV or 
■q'nV . Therefore, terms like 77/77' — > irV followed by V — > ir'y do not contribute to the decays under 
examination. 
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Changing from r)$, rjo to 77, 77' is performed by a rotation involving the (wave-function) 
mixing angle : 



V 




cos 9p 


— sin Op 






v' _ 




sin 9p 


COS #p 







There is, obviously, no loss of generality in introducing this definition and, thus, the 
physical parameter Op which has to be fixed or fitted. 

For all expressions in this Section, the fields which occur are the renormalized ones. It 
should already be noted that all basic Lagrangian pieces involved in the considered 77/77' 
decays do not depend at leading order on the SU(3) breaking mechanism (the parameter 
z = [fx/ 'fir} 2 aheady met), at least in the limit of Isospin Symmetry where we stand. 
However SU(3) symmetry breaking is hidden inside Op (see Eq. ljlOJl ). 



6.2 Amplitudes and Chiral Limit Properties 

With the information given just above, it is an easy task to compute the amplitudes 
for the 77/77' decays considered. One finds : 



T(X — > 7T + 7T 7) = Cx- 



3m 2 n 
1+ p 



D P {s) 



s^e^Kp+Pp , (18) 



using obvious notations, with X = 77,77', m 2 p = ag 2 f 2 and s being the dipion invariant- 
mass. The cx are given by : 

Cf, = [cos Op — xy/2 sin Op] , 

(19) 



sin cos p] 



Finally D p (s) is the inverse p propagator which can be written 

D p (s) = s-m 2 p -U pp (s) , (20) 

in terms of the already defined (KSFR) p mass and of the p self-mass. The occurence of 
Hpp(s) permits to move the p pole off from the physical region (flUEI]- For the present 
purpose, one has only to stress that the p self-mass can be chosen rigorously such that 
n pp (0) = as expected from current conservation [63J. 

Going to the chiral limit, Eq. (|18j) is nothing but the second Eq. (j2J) and has for 
coefficients : 

E 'i(°) =~ A l e 2 J cos0p-xV2smOp] 

(21) 

E> A°) =~ * e 2 \sm0p + xV2cos0 P } . 

These relations clearly meet expectations of Current Algebra ; they could not be de- 
rived exactly if n pp (0) 7^ 0. The single symmetry breaking parameter occuring manifestly 
is x which essentially measures departures from Nonet Symmetry. 
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It should be noted that there is no obvious connection between the way symmetry 
breaking occurs for the box anomaly (Eq. (|21|)) and for the triangle anomaly (see Eqs. 
dHJ) within the broken HLS Model. It is worth recalling once more, that the symmetry 
breaking scheme defined in Section 0] was shown (32] to match perfectly all expectations 
of EChPT collected in (21131120] and that no further piece has been added in order to 
derive Eqs. (|21j) . 

Therefore, the second assumption which underlies the traditional way of breaking 
symmetries for this set of equations (see the discussion after Eqs. ©) is also not met by 
the BKY breaking scheme [39 | PT) | l32j. 

Eqs. (|21|) are also interesting for other aspects : In order to recover the values expected 
for both E' x (0), the VMD (i.e. resonant) contribution happens to be 3 times larger than 
the contact term (CT) contribution and carries an opposite sign ; this was expected when 
building the anomalous HLS Lagrangian Eq. (jSJ). 

Eqs. (118)1 also show that fitting the rj/rf invariant-mass spectra with a constant term 
interfering with a resonant term is indeed legitimate. However, it is also clear that the 
value of this constant is not the value of the full amplitude at origin and thus carries only 
a part of the box anomaly value. 



The triangle and box anomaly expressions in the broken HLS model are summarized 
by Eqs. (jllj) and (J2"T|) . We know from previous works ^HIHSj that experimental data 
support this in the triangle anomaly sector (AVP and 77/77' — > 77) ; the real issue is to 
test its validity in processes where box anomalies are expected to occur. 



6.3 The 77r + 7r tt° Amplitude 



Even if outside the main stream of this paper, it is interesting to give the amplitude 
for the 77r + 7T~7r° anomalous coupling. Using Eqs. (0) and $11}) above, together with the 
piece analogous to Eq. (JTSjl for p ± which can be found in [10], one gets : 



Tin 



7T 7T 



.-_o 



vr u ) 



A(s,t,u) e^e^tPZPp 



[22) 



with 



A(s, t, u) 



1 + 



m. 



D P (s, 



+ 



D P (t) 



+ 



m. 



D p (u) 



(23) 



where m 2 p 



ag 2 f% and D p is given by Eq. (|20|) with the appropriate permutation of the 
argument s to t and u. In this case, the symmetry breaking mechanism we have defined 
has no influence. The momentum dependence of this expression differs from known ones 
(recalled by Eqs. (39) and (40) in 8J) by its taking into account the p self-energy (see 
Eq. pD)l). It gives the expected value (— ie/An 2 fl) ats = t = u = and is worth 
to be checked on forthcoming experimental data [HI]. It might also be extracted from 
e + e~ — > 7r + 7r~7T° annihilation with data covering the low invariant-mass region. 
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6.4 Invariant— Mass Spectra And The Box Anomaly 



^From the amplitude in Eq (fTHjl. one derives the decay width : 

2 

kU , (24) 

for each of the 77 and rj' mesons. We have defined fc 7 = (m x — s)/2mx, the photon 
momentum in the X rest frame and q n = ^ s — 4m^./2, the pion momentum in the dipion 
rest frame. 

The contact term generated by the second piece in Eq. (j5J), is represented in Eq. (I24j) 
by the number 1 inside the modulus squared. On the other hand, as the normalization of 
the VMD contribution can be fixed [4~H l4*2*] at the appropriate value by only normalizing 
to the P77 amplitude in Eq. (|T3|) . checking the effect of this contact term by switching 
on/off this "1" in Eq. (f2"4*|) is indeed meaningfull. In this way, one can address the 
experimental relevance of Eq. (JHJ). 

Eq. ()24|) is interesting in many regards : 

• The shape of the invariant-mass spectra depends on the 77/77' meson properties 
only through a kinematical factor (k^). Therefore, the shape of the invariant-mass 
spectra does not carry any manifest information on the box anomaly constants cx, 

• The lineshape of the invariant-mass spectra in 77/7/ decays depends only on p meson 
properties. However, the way this dependence occurs in 77/77' decays is different from 
the one in the pion form factor |4H] . as the dressing of the p — 7 transition amplitude 
n p7 (s) plays no role in the 77/77' decays, 

• All information on the value of c v and c v > is carried by the partial width itself. It 
can be algebraically derived if D p (s) is known reliably from another source. 

In order to perform a search for the box anomaly, one needs a function D p (s) ac- 
curately determined between the two-pion threshold and the mass. In the physical 
region involved in 77/77' decays, all coupled channels allowed by the HLS model contribute 
at one-loop order |4~3~] . However, except for 7r + 7r~, each provides only logarithms, beside 
their influence on the subtraction polynomial hidden inside the p° self-mass U pp (s). This 
is their major effect, and thus neglecting these loops while still considering a subtraction 
polynomial of the appropriate degree is certainly motivated. 

Reduced to only its coupling to 7r + 7r~ (with or without accounting for kaon pairs), the 
p propagator used here contributes to providing a fairly accurate numerical determination 
of the pion form factor both in modulus jlSllH] and in phase [13] up to the <fi mass. 

7 And, to some extent, except also for the lott q channel in the rf decay ; however, the neglected effect 
is concentrated in the region above 917 MeV, very close to the phase space boundary for 77' decay and 
far beyond in the rj decay. 



dT(X — > 7T + 7T 7) C 



X 







d^fs 



36 [27TA] 



1 + 



3m; 
DM 



14 



Therefore, for the purpose of studying the box anomaly, there no point in going beyond 
contributions from only the non-anomalous HLS Lagrangian. In this case, the p self- 
energy can be written : 

n«»(«) = 4[^) + ^M«)] » (25) 

where g pn7T = ag/2, while z has been already defined. We have denoted &k{s) and £k(s) 
the pion and kaon loops amputated from their couplings to external legs (we neglect the 
mass difference between K ± and K°) ; these are given in closed form in [12] , 

These loops should be subtracted at least twice in order to make convergent the 
Dispersion Integrals which define them as analytic functions ; this gives rise to a first 
degree polynomial P p (s) with arbitrary coefficients to be determined by imposing explicit 
conditions or by fit. However, as noted just above, anomalous loops force to perform, 
at least, three subtractions [I2HM], which modifies the arbitrary polynomial P p (s) to (at 
least) second degree. It is the reason why P p (s) will be assumed of second degree, even if 
one limits oneself to pion and kaon loops. This does not increase our parameter freedom, 
as will be seen shortly. 

6.5 External Numerical Information 

^From what seen above, the condition P p (0) = on the subtraction polynomial is 
certainly desirable, otherwise the Current Algebra expectations could not be derived 
exactly ; additionally, this condition ensures masslessness of the photon after dressing by 
pion and kaon loops. It thus remains 2 subtraction constants to be determined or chosen ; 
we shall fix them from the fit performed [E] on the pion form factor from threshold to 
the mass. Denoting U pp (s) the p self-energy subtracted three times 02], we have : 

II pp (s) = TT pp (s) + ei s + e 2 s 2 . (26) 

On the other hand, we can also fix the HLS parameters a, g (and thus m p ) and x 
to their values fitted in radiative and leptonic decays 8 |5TJ|12]- Knowing z and x, one 
can derive the value for 9p from Eq. (jlOj). The values for e.\ and e 2 are fixed from an 
appropriate fit to the pion form factor [33] , where the parameters a, g and x are fixed 
consistently with AVP and (uj/<j))e + e~ modes. 

As we restrict the p coupling to only the 7r + 7r~ and KK channels, these values for e% 
and e 2 are certainly correlated with the chosen values for a and g ; they are not affected 
numerically by the value for x. 

The values for these parameters are gathered in Table H As these parameters are 
supposed universal in the realm of the HLS Model, one can fix their values from fit to 

8 Actually, the values for g and x are determined almost solely by the AVP radiative decays ; the 
value for a is a consequence of these on the V — > e + e~ decay modes, but mostly the u> and <f> leptonic 
decays. 
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Parameter 



Value 



a 2.51 ±0.03 

g 5.65 ±0.02 

x 0.90 ±0.02 

z = [f K /Uf 1-51 ± 0.02 

~ 0.222 ±0.011 

_e2 -1.203 ±0.017 GeV" 2 

Table 1: Parameter values for a, g, x ,z fixed from a global fit to VP'j and V —>■ e + e~ 
decay modes |Hll32lEn] ; &i and e<i are fixed from a fit [13] to the pion form factor 
including only tttt and KK as channels coupling to p. 



data independent of the 77/77' — > 7r + 7r~7 decay modes. It is indeed the case for the VP^y 
or the Ve + e~ decay modes and for the pion form factor. As commented on above, these 
fit values correspond to a very good fit quality for the corresponding data. For instance, 
they allow to predict the two-photon decay widths as recalled in Table El 

Parameter PDG 2002 Prediction Significance (n a) 

^ 77 (keV) 0.46 ±0.04 0.46 ± 0.03 0.0 a 

7/->77(keV) 4.29 ±0.15 4.41 ± 0.23 0.4 a 

Table 2: Partial widths for 77/77' — > 77 as predicted using Eq. (jllj) with parameter values 
as coming from a global (HLS) fit to only VP7 decay modes of light mesons [3*2"] . 



Choosing the p propagator as it comes out of the HLS fit to the pion form factor [43J is 
also fairly legitimate, as this p propagator should also be valid anywhere within the HLS 
framework. As stated above, we consider for clarity only the case where the only open 
channels are nir and KK. We have, nevertheless, checked that changing to various open 
channel subsets coupling to the p meson (as done in jlHj for the pion form factor), with 
correspondingly changing e± and €2 to their fit values, does not produce any significant 
modification to the results presented below. 

To summarize, self-consistency implies that we can fix all parameters and functions 
from their most reliable fit values and expressions, provided the data set is independent of 
the 77/77' decays considered here. This independent data sample covers VPj and Ve + e~ 
couplings |4*T | l4*2" l l3*2*] and the pion form factor [13! • ^ then follows that all information 
related with the box anomalies can be predicted without any parameter freedom. 
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7 Predictions For tj/t/ —> n + n 7 Decays 

In this Section, we examine the predictions derived for the 77/7/ — > 7r + 7r _ 7 decay 
modes for the partial widths and dipion invariant-mass distributions. 

7.1 Partial Widths, Experimental Values and Predictions 

Using Eq. and numerical (and functional) information given in the previous 

Subsection, it is easy to check that we can write : 

A x c\ , (27) 

where : 

k*q%dyfi , X = 7 l ,rf . (28) 



2 

1 a 

A, 



r(x 


— > 7T + 7T 7) 










J 2m.Tr 


D p (s) 



36 [27r/ ff ] 

This integration can be done by Monte Carlo techniques and gives : 

A v = 38.25 ± 1.07 , eV A v , = 42.16 ± 3.00 keV . (29) 

For further concern, one should note that these integrals are not affected by the value 
for x. Using the parameter values given in Tabled Eqs. and (|1()|1. one can compute 
the partial widths and get the results collected in Table El We have performed the 
computation by switching on/off the contact term contribution 9 . We stress again that 
all results presented in this Section do not depend on any free parameter and thus are 
predictions relying on the rest of the HLS phenomenology. 



decay 


PDG 2002 


Prediction with 
Box Anomaly 


Prediction without 
box Anomaly 


i] — > 7r + 7r~7 (eV) 


55 ±5 


56.3 ± 1.7 


100.9 ±2.8 


Significance (n a) 




0.25 a 


8 a 


t]' — > 7r + 7r~7 (keV) 


60 ±5 


48.9 ±3.9 


57.5 ±4.0 


Significance {n a) 




1.75 a 


0.39 a 



Table 3: 77/7/ Partial widths as predicted by the HLS Model when switching on/off the 
box anomaly contribution. Significance is computed using an error obtained by adding 
in quadrature the experimental error and the relevant model error computed by Monte 
Carlo sampling (using information in Table P) . 



9 When switching off the contact term in Eq. J33J), the numbers in Eq. (|29|l become A v — 57.51 ±4.01 
eV and A v i = 49.60 ± 2.98 keV, which is already conclusive. 
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^From Table El one clearly sees that the predicted partial width for rf is not really 
sensitive to the presence of the contact term. This can be well understood as, indeed, 
the value for A v i is sharply dominated by the p peak contribution provided by the VVP 
Lagrangian term and the magnitude of the contact term is comparatively small. 

In contrast, the predicted partial width for r\ is much more sensitive to the contact term 
because this contribution has only to compete with the low mass tail of the p distribution ; 
the bulk of the resonance contribution is indeed sharply suppressed because the available 
phase space is small and located far outside the p peak. 

Therefore, one can already conclude from Table El that the T]/r]' partial widths values 
provide a strong evidence in favor of the box anomaly. Unexpectly, this evidence is 
provided by the rj partial width alone. Additionally, the values predicted for the box 
anomaly constants c v ~ 1.21 and c v > ~ 1.07 from the rest of the HLS phenomenology 
fits nicely the rj/rj' partial widths, which means, for instance, consistency with having 
6 P = 10.30° ±0.20°. 

Together with the results predicted for two-photon decay widths of the rj/rj' system, 
this also gives a strong support to the extended BKY breaking scheme summarized in 
Section |U and to Eqs. (fTTj) and (}2*Tj) for the amplitude expressions at the chiral limit. 

7.2 Invariant— Mass Spectra With/ Without The Contact Term 

The shape of the dipion invariant-mass distributions are given in Figure ^ top for 
rj', mid for rj. These are proportional to yields (up to acceptance/efficiency effects). The 
distributions are displayed with having switched on/off the contact term ; in these two 
figures, the relative magnitude of the twin distributions is respected. 

Looking at the 77 distributions, one clearly understands the width results given in 
Table EH as the integrals corresponding to box anomaly on/off are clearly very different 
(actually by a factor of about 2). 

In the case of the 77 meson, lineshape differences between the case when the contact 
term is activated and when it is dropped out are tiny as illustrated by Fig. Q bottom. 
In this figure, one displays the distribution obtained by removing the contact term and 
the one derived by activating it, after rescaling it by ~ 1.8. 

The lineshape, in the case of the rj', shows that the peak location when accounting 
for the contact term is slightly higher mass (6-7-8 MeV) compared to the case when this 
(CT) contribution is cancelled out. However, the main effect is that yields below the p 
peak location are somewhat suppressed because of the contact term. 

In this sort of situation, if one performs a fit of an rf spectrum affected by CT with 
only a resonance contribution and lets free the resonance parameters, the shape will be 
distorted. Indeed, in order to average reasonably the rising wing of the p distribution, 
the peak has to be shifted to higher mass and therefore the observed mass must be larger. 
This is a mechanical effect connected with the minimization of a \ 2 f° r an y appropriate 
function of one variable. We come back to this point when comparing with experimental 
data. 
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8 Experimental Data on rj/rf —> n + n 7 Decays 

There are several sets of data available for the dipion mass spectrum of the r( meson. 
Most of them have been published only as figures OUHlEHIISIinilHl • For these, how- 
ever, it happens that the information given in the body of the articles provides enough 
information in order to recover the yields and derive the acceptance/efficiency function ; 
the redundency of the information is fortunately such that consistency checks and cross 
checks can be performed which validate the outcome of the procedure. This is described 
in details in Section 4 of [33] together with the peculiarities of each of these data sets. 

Other spectra were available directly to us [TH] or as PhD theses (12112011211122] pub- 
lished only as preprints ; Here also the relevant information was either directly available 
or could be reconstructed accurately, as for the references quoted above. One should note 
that the data of (21] supersed the ARGUS results published in [T4] . 

These former data samples carry widely spread statistics ; 474 events for the oldest 
data set [12], 130 events from TASSO [13], 795 events from ARGUS [H] updated three 
years later to 2626 [21], 321 events in the TPC-77 sample [15], 195 events in the PLUTO 
data set PU, 586 in the CELLO data set [22], 401 for the data set of WA76 collected 
using the Omega Prime Spectrometer at the CERN SPS ^H] and, finally, 2491 (after 
acceptance corrections) for the experiment performed at Serpukhov using the Lepton F 
facility [H]. 

The method used to extract the dipion invariant-mass spectra from data is of special 
concern. These were derived from the data samples just listed in the following way : for 
each bin of dipion invariant-mass, one plots the 7r + 7r~7 invariant-mass spectrum and fits 
with a gaussian (plus a polynomial background) the number of 7/ it contains. In this way, 
one get rid to a large extent of the precise background 10 parametrization, as the signal 
is a narrow gaussian peak dominated by the experimental resolution. 

Performing this way, spectra appear without any background and the influence of this 
in the data sets only reflects in the magnitude of the errors on the yields per bin. It should 
be stressed that this extraction method is obviously independent of any assumption on 
the lineshape of the underlying p invariant-mass distribution. 

For the data samples of [TTlIIHlUn] , the acceptance/efficiency function was directly 
known, but without information on its uncertainties. As the spectra of [TS1E3 carry 
small statistics, statistical errors are dominant and errors on the acceptance/efficiency 
function can be neglected. For the data sample of ^7j, the acceptance/efficiency function 
is provided as a curve (see their Fig. 3) ; as no information is reported about uncertainties 
affecting this function, these cannot be accounted for when folding in this function with 
any model distribution. 

For the other data samples reviewed above, uncertainties on the acceptance/efficiency 
functions are also unknown, as these can only be derived by unfolding it from the fitting 

10 We have to make assumptions on the background shape across some small 7r + 7r~7 mass interval while 
the signal is a narrow gaussian (typically 20 to 30 MeV for its standard deviation). This is certainly 
much safer than assumptions on the background shape over a 1 GeV invariant-mass interval with on 
top of this a signal as broad as a p distribution. 
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distribution. This was always provided as a product of a well defined model function for 
the decay with this acceptance/efficiency function. This is also of little importance for 
all data sets dominated by statistiscal errors, but it also affects the large statistic sample 
of ARGUS we shall examine [2*Tj . 

Neglecting this source of uncertainties when computing model errors certainly biases 
X 2 estimates towards larger values (and smaller probabilities). However, it should not 
spoil qualitatively model descriptions. 

The sample of Markll [16] is also significant (~ 1200 events), however, the mass 
spectrum derived from this has been obtained in a different way : Selecting the events in 
some mass interval around the rf mass in the global 7r + 7r~7 invariant-mass spectrum, the 
corresponding events are plotted in bins of 7r + 7r~ invariant-mass. This spectrum is then 
described as a superposition of a p mass distribution plus some background, and a global 
fit to this spectrum provides the signal (p) and background populations inside each bin. 
Therefore this method assumes an accurate knowledge of all phenomena contributing to 
the background (and of its parametrization) ; it also relies on the way the p lineshape is 
parametrized. This is also, basically, the method used to study the rj' mass distribution 
performed by the L3 Collaboration [21] on a sample of 2123 ± 53 events ; this will be 
specifically discussed at the appropriate place below, as it is the latest published data 
sample. 

Some other papers published spectra without background subtraction (namely [02] 
and £2,Q] which carry actually the same data). In order to use these, one would have to 
model the background without any motivated knowledge of the data set and detector 
properties 11 ; therefore, this Markll spectrum will not be examined here. This lack 
of background subtraction is also the reason why the spectrum published by the L3 
Collaboration is also skept. 

Finally, the most reliable spectrum for the rj' decay is the one collected by the Crystal 
Barrel Collaboration [23] which is also, by far, the largest data sample (7392 events). 
This spectrum has been constructed using the method described at the beginning of this 
Section ; therefore, it is independent of any assumption on the underlying p invariant- 
mass distribution and, thus, is certainly free from any prejudice or bias. Additionally, 
this data sample is certainly the most secure to be used as uncertainties on acceptance 
and efficiencies are already included into yield errors and, thus, model comparison can 
be performed directly and reliably. 

The corresponding spectra for the 7] decay, have been derived from up — > rjn data 
collected long ago by two experiments [HHHT] ; the published data are already background 
subtracted. They both carry large statistics (7250 for jTU| and 18150 for [H]). The 
relevant results have been published only as figures. Yields per bin can be read off from 
these figures without any difficulty together with their errors. 

In jTUj, the acceptance/efficiency function is provided directly and also folded in with 
a well defined model function. It is given superimposed to the case when the analysis 

11 Indeed, beside extracting the yields, one needs to estimate the acceptance/efficiency function which 
might well be different for signal and background events. 
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is performed with a simple gauge-invariant phase-space matrix element and with a p 
dominant one. The two corresponding acceptance/efficiency functions are conflicting, 
essentially in the region fc 7 = 90 H- 110 MeV. However, this seems to reflect their depen- 
dence upon angular distributions. It is therefore worth using the functional information 
associated with the p dominant matrix element. 

For the data of Layter et al. llj, the acceptance/efficiency function is not shown 
and should be unfolded from the theoretical p (and phase-space) distribution(s). This 
information can be extracted with some reliability ; in contrast with |10j . this yields 
a function extracted from the p distribution very close from those extracted from the 
simple phase space distribution. Actually, this data set should be considered with some 
care as extracting the acceptance/efficiency function can only be performed by making 
some assumption on the p mass actually used in this paper ^1] . We have conservatively 
assumed that Layter et al. ^1] used the same p mass as Gormley et al. [IjQJ, namely 
m p = 765 MeV ; this assumption is crucial and cannot be ascertained. This makes more 
secure information derived from the Gormley spectrum. 

Finally, for both t] spectra, it is impossible to restore the accuracy on the accep- 
tance/efficiency function. These will be considered negligible in the present study. 

9 Experimental Data Versus Predictions For 77/77' 
Spectra 

As seen in Subsection 16.41 the HLS Model provides definite spectra for both rj/r]' 
invariant-mass distributions. These depend on parameters which can be fixed indepen- 
dently of the r]/r]' — > 7r + 7r _ 7 spectra, like a, g, x (see Subsection 16. 5|) . and of the p 
propagator which is fitted elsewhere 03] with parameters values as determined in these 
fits (see Eq. (|2T)j) and Table [T]). The model fairly predicts the absolute magnitude (the 
integral) of each spectrum as illustrated in Table H3 In this Section, we focus on com- 
paring the predicted lineshapes derived from the model Eq. (|2~4"|) with the data listed in 
Section |S1 

As all data considered are binned, we have integrated the predicted function Eq. 
(J2U (or Eq. (HEI)) over the bin size and normalized this function to the integral of the 
experimental distribution. When relevant (i.e. all spectra except for the one of Crystal 
Barrel j2Hl)> the model function was folded in with the acceptance/efficiency function 
derived for each of the above data samples. 

The results are displayed in Figs. (j2J) to (JU). One should note that all rj' spectra are 
given as functions of the dipion invariant-mass, while 7] spectra are given as functions of 
the photon momentum in the rj rest frame. 

In these figures, together with the specific experimental spectrum, we show the pre- 
dicted curve (computed just as defined just above) when keeping the contact term (full 
curve) and the one derived by dropping out this contribution (dashed curve). These 
two cases will be referred to as resp. CT and NCT. In these figures, we give the \ 2 
corresponding to these two solutions under the form x 2 (CT)/x 2 (NCT). The number of 
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degrees of freedom can be easily read off from the spectra as this is exactly the number 
of bins of the experimental histogram. 

The curves shown have been computed at the central values of the parameters as given 
in Table ["[] The x 2 's have been computed by folding in the experimental error in each bin 
with a model error also computed bin per bin. These model errors have been computed 
by sampling the parameters around their central values with standard deviations given 
by their quoted errors (see Table Q). Except for Crystal Barrel [23], where it is irrelevant, 
uncertainties on the acceptance/efficiency functions are not (cannot be) accounted for. 
The curves shown are actually histograms which have been smoothed automatically by 
the hbook/paw package. 

Examination of Figs. (J2J) to (""Q) is quite interesting. First of all, the spectrum from 
TPC-77 is clearly the single one far away from predictions, indicating that something 
was not well controlled when extracting it from data. All others match well, or quite 
well, the predictions ; this clearly gives support to the model developped in the previous 
Sections and to the relevance of the parameters given in Table [*[J 

In terms of probabilities (reflected by the y^jdof values given in the figures), the 
oldest data set of [12] gives comparable probabilities to either of the CT/NCT assump- 
tions, while maybe slightly prefering the NCT assumption. CT/NCT descriptions are 
practically equivalent for the ARGUS [21] and WA76 [TS| spectra, while nevertheless 
slightly favoring the CT assumption. 

The relatively low statistics spectra provided by TASSO [13] (x 2 ratio of 0.6 in favor 
of the CT assumption), CELLO [22] (0.7) and PLUTO PH (0.7) somewhat prefer the 
CT assumption. 

Finally, the two largest statistics experiments Lepton F [17] and Crystal Barrel [2*3] 
sharply favor both the CT assumption against NCT ; the \ 2 distance is indeed better 
by more than a factor of 2 . 

To be more precise the Lepton F spectrum 12 gives a 3% probability to the CT as- 
sumption and a 2 10~ 8 % probability to the NCT assumption. These relatively low 
probabilities should be related with the lack of information on the acceptance/efficiency 
function which affects in a same manner both solutions. Accounting for the corresponding 
errors would certainly increase both probabilities but hardly switch their ordering. 

The corresponding probabilities for the Crystal Barrel data set are resp. 57.8% and 0.1 
% ; these values are certainly realistic as the model errors are reasonably well accounted 
for. 

Among these two data sets which could be used as reference, the Crystal Barrel data 
(available in a directly usable form [23]) should clearly be prefered, as systematics are 
better controlled all along the invariant-mass range. It also carries, by far, the largest 
statistics. 

12 The 8 lowest mass points of this spectrum contribute severely to the \ 2 f° r both (CT/NCT) 
assumptions. On the other hand, the sharp drop in acceptance |17j at large might have been 

difficult to estimate reliably. Qualitatively, however, the clear preference of this distribution for the CT 
assumption is obvious. 
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^From x 2 values, the shape of the n spectrum from Layter et al. jTT] seems in better 
agreement with the NCT assumption (77% probability) than with the CT assumption 
(3% probability). 

In contrast, the description of the r\ spectrum from Gormley et al. ^Uj is simply 
perfect and corresponds to resp. 97.7 % probability for the CT assumption and to 58.2 
% probability for the NCT assumption ; this reflects better the remark following from 
Fig. (bottom) that these lineshapes are very close together. The \ 2 values, however, 
indicate that the geometrical distance (~ y/x*) of this spectrum to the CT solution is 
significantly smaller than those to NCT. 

Fig. El gives the same information as in Fig. but enlarged and binned. Here one 
sees that a third of the x 2 f° r the Layter spectrum ^T] comes from only the bin covering 
the momentum interval 60 -j- 80 MeV/c. Compared to the same result for the Gormley 
spectrum [TO], the Layter spectrum looks a little bit skewed. It is, however, impossible to 
decide whether this comes from systematics affecting the acceptence/efficiency function 
as this (skewed) shape happens to match nicely the NCT assumption 13 . 

However, the rj — > 7r + 7r~7 partial width alone certainly a more secure information, 
and the Gormley spectrum undoubtfully favor the CT assumption against the NCT one. 
These two aspects have to be balanced in a global fit accounting for lineshapes and partial 
widths. 

10 A Global Fit to 77/77' Spectra and Widths 

We have compared the data (lineshapes and partial widths) with the predictions of 
our model fed with numerical and functional information coming from the rest of the 
phenomenology accessible to the HLS framework, without any parameter freedom. The 
results obtained in Section and El considered together indicate that the model is valid 
and favors the the contact term as a physically motivated contribution to decay processes. 
We remind that this contact term is not a free parameter, as widely discussed above. 

In view of this, it looks worth performing a simultaneous fit of the n data sets with 
some accurate n' spectrum ; for reasons explained above, it is certainly worth choosing 
the Crystal Barrel spectrum. As a clear conclusion should take into account all aspects 
of the available experimental information, partial widths have been fed into the \ 2 to be 
minimized. 

In order to perform this fit one needs to release some of the parameters fixed as 
in Table ^ ; as the main information for the present purpose is the peak location, it 
looks worth releasing the parameters named e\ and €2 which mostly influence the p 
peak location. Comparing the values returned from this fit to the corresponding values 
originally extracted from fitting the pion form factor could contribute to clarify the 

13 This skewness might have been magnified unwillingly by the choice of the m p value we performed in 
order to extract the acceptance efficiency /function for the Layter spectrum. Any underestimation of this 
input nip value contributes to the skewness of this distribution. We are responsible for this uncertainty, 
but we did not find an unbiased way out. 
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conclusion, as one can consider the p propagator as a universal function, as valid for 
F n (s) as for the the rj/rf spectra. 

Indeed, we know that, in the pion form factor, the subtraction polynomials of U pp (s) 
and il p7 (s) are somewhat competing and that some (small) correlation among the cor- 
responding polynomials exists (43 j ; it is therefore motivated to attempt freeing e\ and 
€2 as these correlations could have spoiled their central values by some (certainly) small 
amount. However, the parameter values returned from fit must not be inconsistent with 
their partners derived from fit to the pion form factor only 14 . 

On the other hand, the other parameter values in Table Q describing fairly well the 
full set of V — > P7 and P — > decays would hardly accomodate a significant change 
of their values without failing to fit the VP'-f processes. 

In order to avoid too much correlations which can hide clarity in the conclusions, 
we shall test separately the CT and the NCT assumptions. Indeed, as the HLS Model 
predicts the magnitude of the constant contact term (if any), it seems enough to check 
its precise relevance and no attempt will be made to fit its value. Finally, for the present 
exercise, we neglect model errors 15 ; this mechanically makes the x 2 's slightly more pes- 
simistic than they really are. 

We could have chosen to perform a simultaneous fit of the Crystal Barrel rf data 
set [23| together with both rj data sets 11. 10J simultaneously. One could indeed imagine 
that the systematics could compensate. We have, nevertheless, prefered performing the 
fits separately for the Crystal Barrel 7/ spectrum together with each of these 77 spectra 
in isolation. Using both 77 spectra certainly leads to intermediate fit qualities. 

Additionally, before letting e% and €2 vary, we have performed the "0 parameter fit" 
in order to get the x 2 ' s an d probabilities when using directly the parameter values as 
given in Table ^ In this way, we know the starting quality of the global description of 
these decay modes induced by the rest of the HLS phenomenology ; we can also estimate 
what is gained by letting some parameters to vary. 

When using the data set of Layter et al. [TT] . while accounting for the 77/7/ contact 
terms at the expected level (assumption CT), one clearly sees from Tabled that the p 
lineshape parameters e\ and €2 do not move farther than 2 a from the values found when 
fitting the pion form factor 03] (the present cr's are, however, much larger than found 
in fits to F n (s) f43 ; ). The gain in \ 2 g°t by releasing these parameters is modest (2.7) 
and the central values for the partial widths get a little changed. This confirms that the 
parameter values in Table [T] giving x 2 /dof = 37.88/35 (34% probability) are already close 
to optimum ; leaving them free, essentially improves the 77 spectrum lineshape slightly, 
but at the expense of slightly degrading the central values of the partial widths. 

14 This actually means that a further test could be a simultaneous fit, within a consistent framework, of 
the pion form factor and of the relevant rj/rj' decay information. One does not expect neither a surprise 
nor hard difficulties from such an attempt ; the present work indicates that this should not provide more 
insight than a global probability. 

15 These are certainly present as the uncertainties on a, g and x contribute to model errors, even when 
releasing any constraint on e\ and e-i. 
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When dropping out the contact term contributions (assumption NCT), the p line- 
shape parameters e% and €2 change significantly with respect to their starting point, with 
much larger errors than originally. Even then, the fitted values for e\ and €2 move by 
more than 6 (new) a from expectations ; therefore, these fitted values can be considered 
inconsistent with their values fitted in the pion form factor. Additionally, even if the gain 
is large (x 2 /dof is improved from 140.65/35 to 60.88/33), it is not sufficient to push the 
probability (0.2%) to a reasonable value. Therefore, the peculiar uncertainties affecting 
this spectrum does not prevent to reach a clear global conclusion. 

When using the data set of Gormley et at. [TU] together with the rf data of Crystal 
Barrel, the picture is unchanged, but looks much clearer. Fixing the parameters to their 
values in Table 1, gives already a remarkable fit quality (probability 85 %), when contact 
terms are accounted for. In this case, letting free e\ and e-ii the x 2 improves a little (0.57 
unit), but the fit probability degrades to 80 %, because of the smaller number of degrees 
of freedom. The p parameters move by about 1 (new) a from expectations and thus stay 
consistent with the p lineshape determined when fitting the pion form factor ^Hj (the 
region of the minimum \ 2 seems flat). 

When removing the contact term from our expressions (NCT assumption), the start- 
ing values of the fit parameters provide a comparatively poor decription (1.6 % prob- 
ability) and the 77 partial width is far from expectations [5| (see also Table 0}. Now, 
releasing e\ and €2 improves significantly the description, as we reach a 27 % probability 
after fit, with reasonable central values for both partial widths. The price to be paid for 
this configuration is that the parameters e\ and €2 change by more than 3 (new) a from 
expectations. Therefore, the NCT assumption returns a p lineshape inconsistent with 
fits to the pion form factor. 

^From the previous Sections, we already knew that the CT assumption is certainly 
favored in a global account of both shape and partial width for both the rj and rf meson 
simultaneously. We also knew that the NCT assumption was disfavored under the same 
conditions. 

What we have learnt in this Section is that, in order to accomodate the description of 
all aspects of the 77/7/ information, the NCT assumption gives up being consistent with 
the p lineshape as found by fitting the pion form factor • 

Therefore, we conclude that experimental data do provide a fair evidence in favor of 
the box anomaly phenomenon at the expected level ; additionally, the sharing observed 
between resonant and contact term contributions (—3 : 1) is well predicted by the FKTUY 
assumption [3U] leading to the Lagrangian in Eq. (jSJ). 

In Fig. we show the description of the Crystal Barrel spectrum [23] using Eq. 
(I24|) (or Eq. (J28)0 with the contact term considered and removed. In order to get this 
we performed fits leaving free e\ and €2- 

When accounting for the vf contact term, the p peak location is found in the bin 
covering the mass region from 725 to 750 MeV. When dropping it out, the (fit) mechanism 
described in Subsection 17.21 makes the p peak shifting to the next bin which covers the 
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mass interval from 750 to 775 MeV. This trend was already observed by f7H2H] and also 
by most Collaborations who have performed extraction of the rj spectrum canonically ; 
sometime too much jT^I- 

A real shift exists and is small (see Section 1772}) . It is artificially increased by the fit 
procedure in order to get a better account of the low mass tail of the rf invariant-mass 
spectrum. However, this artificial large mass shift is indeed the signal of the box anomaly. 

It is claimed in [21] that the L3 Collaboration does not observe a p peak shift. Several 
reasons can be invoked. First, as remarked above, [21] did not perform the rf spectrum 
extraction canonically and, therefore, any conclusion about the underlying p lineshape in 
the rf decay becomes a delicate matter. 

Among other reasons the most likely is their fitting of the p mass and width. The 
p (Breit-Wigner) mass is thus found at values (766 ± 2 MeV) normally obtained in 
only processes where the dynamics is not really well under control (hadropro duct ion or 
photoproduction) 

If, for a moment, the L3 result were considered as reference in order to detect a p 
peak shift, one might have instead to consider that a shift occurs in e + e~ annihilations or 
t decays, as these yield rather larger p (Breit-Wigner) masses (~ 775 MeV) jH]- Under 
these conditions, it is difficult to draw any conclusion from [21] about the existence (or 
absence) of a p peak location shift in the rf — > 7r + 7r~7 decay. 

11 Fits To The Four Anomaly Equations 

/,From now on, we make the assumption that the correct set of equations defining the 
anomalous amplitudes at the chiral limit are given by Eqs. and (j21J) and no longer 
by Eqs. Q. These have been derived using the approximate field transformation Eq.fjHJ). 
We also examine, for completeness, the case when the exact field transformation is used ; 
the way to modify our anomaly equations to go from one case to the other is given in the 
Appendix. 

In both cases, these equations actually depend on only one parameter (resp. x or A) ; 
this can legitimately look like a severe constraint. 

11.1 Fit Results With Approximate Field Transformation 

These equations depend only on f w , z = [fx/ fw} 2 and on x. In the present framework, 
9p is no longer an independent parameter as it can be algebraically derived from Eq. (jl(J|) . 

One can consider legitimate to still fix f n to its experimental value (92.42 MeV) ; this 
is also true for z (see Table Q). Therefore, our set of anomaly equations depends on only 
one parameter x, we choosed previously to fix from fit results to radiative decays [4*Tjl32| . 
Releasing the constraint Eq. (fTUj) would only add a comfortable (and useless) parameter 
freedom to the fits presented just below. 
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Therefore, one considers here Eqs. and (}2*T]) by themselves and attempt to fit 

them as a constrained system of 4 equations with only one unknown (x). The results are 
expected to provide consistency with those obtained for the same parameters and physics 
quantities derived elsewhere ^UE21 from fit to VP'-f decay modes. 

The 77 partial widths are related with the amplitudes given in Eqs. by : 

r(X-77) = ^||Gx(0)| 2 (30) 

On the other hand, the partial widths T(X — > 7r + 7r~7) given by Eq. ()27j) where the 
coefficients Ax given by Eq. (J2HJ), depend only on p properties already derived in by 
a fit to the pion form factor. The errors on Ax are taken into account in the fit as they 
are independent of x. 

One has performed a fit of these four partial widths keeping first z fixed and allowing 
x to vary. The results are summarized in Table 03 

It is clear that the fit is fairly successfull and represents the most constrained fit of 
the four partial widths ever proposed. One should remark that the best fit returns a 
value for x perfectly consistent with our previous fits to solely the VP^y radiative decays, 
as can be concluded by comparing to its input value (see Table HJ). The corresponding 
value for 6p is not changed compared to our previous estimates from fit to VPj decay 
modes : 9 P = -10.48° ± 0.18°. 

We do not give the estimates for derived quantities (fo, f$, , 9s) as they practically 
coincide with the values given in [32] and are all in good correspondence with expec- 
tations. Concerning partial widths, three out of four reach a significance much better 
than the 1 a level ; the worst case is Y{rf — > 7r + 7r~7) for which the distance to the 
recommended value 9] is "only" ~ 1.6 a. 

The fit quality yielded (x 2 /d°f — 2.66/3) is such that releasing also z can look like 
an academic exercise. It has nevertheless been performed as some correlation could spoil 
numerically the connection between z and [fx/fir} 2 - 

The fit returned x = 0.908 ± 0.021 and z = 1.488 ± 0.054 with x 2 /dof = 2.62/2, 
practically unchanged, corresponding to a 27 % probability (the worse significance is 
due to having less degrees of freedom). The correlation coefficient is +0.67, and the 
minimization does not spoil the numerical values found elsewhere |4*Tll3*2*] for the same 
parameters. 

Therefore, this leads us to conclude that the VPj decay modes on the one hand, 
and the four standard anomalous r]/i]' decay modes on the other hand, yield information 
fairly consistent with each other. This also means that the anomaly equations we derived 
are consistent and that the approximate field transformation (leading order in breaking 
parameters) on which they rely match well the present level of accuracy of the data. This 
statement will be confirmed directly shortly. 
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11.2 Op versus x 



Eq. (JTU|) corresponds to setting the EChPT decay constant OE] F® to zero. This is 
rigorously expressed in the broken HLS model by : 

tmSp= (WW) = A (31) 



in term current matrix elements and of their expressions [22] • O n the other hand, detailed 
computation yields : 

f = H±£ , £ = f (l-,)x (32) 

Jit OX Jtt O 

As clear from its expression /o keeps the first non-leading contribution in breaking 
parameters (x — 1 + [x — 1]) ; for 6 one has naturally chosen to replace x by 1 and this 
leads to Eq. (jTUjl . However, one may be tempted to keep it and this leads to replace x 
by x 2 in the expression Eq. (|10|) for tan#p ; this is nothing but changing the existing 
term of order 0([z — l][x — 1]) ~ 0.05. 

We have redone the fit just described with this change and yielded a slightly better fit 
quality than the previous one (x 2 fdof = 2.61/3). This fit returns also x = 0.902 ± 0.017 
(~ 0.5 a from its partner in Table El or also a one percent change) and no change at all 
for the partial widths compared to what is displayed in Table ED 

Therefore the sensitivity in describing data is not sharply dependent on non-leading 
contributions in tan#p and using Eq. (JTUJl with x or x 2 gives undistinguishable results, 
while the latter might be prefered. 



11.3 Fit Results With Exact Field Transformation 

In order to check the sensitivity of the model to some other details of the broken HLS 
model, we have also attempted fits using the exact field transformation 32j instead of its 
leading order approximation (see Eq. (JSJ) ; some details and formulae are given in the 
Appendix. 

The main motivation was to figure out the sensitivity of the data to the approximation 
performed on the field transformation. 

In this last series of fits, we have kept z fixed ; therefore the fitting parameters are 
A (the basic Nonet Symmetry breaking parameter, see Eq. (J7J)) and P) the later being 
possibly fixed by the constraint 0q = 0. The fit results and physics quantities of relevance 
(ChPT decay constants, mixing angles and partial widths) are given in Table El 

The first conclusion one can draw from this last Table is that the fit value for 9q departs 
by less than 1 a from zero and the consequences of this on derived physics quantities is 
simply negligible. Stated otherwise, the present data are insensitive to releasing the 
constraint #o — 0. This constraint allows to extract a value for 8g = —18.2° with a 
very small statistical error (~ 0.25°) ; The value found for fo and fg are in the usual 
ballpark and nothing noticeable appears compared to the case when the approximate 
field transform was used |32|. 



28 



The partial widths are still quite consistent with those in Table showing that the 
refinements introduced by the exact transformation have no impact on the extracted 
width information for 77/77' — *> 77 and 77/77' — > 7r + 7r~7. 

One might maybe note that, in all our attempts, we never get a solution with the 
partial width for 77' — > 7r + 7r~7 larger than its recommended value jH] ; so, the observed 
1.6 a departure looks like some small systematic effect. This could be due to having 
neglected some unidentified tiny (higher order) contribution ; this might also indicate 
that the recommended value is slightly overestimated. 

On the other hand, we have also reconsidered the problem of which value for 77 — > 77 
should be prefered among the the recommended value f5] -recently confirmed by a direct 
measurement of this branching fraction [65j-. the 77 measurements and the (single) 
Primakoff effect measurement. This was done already in j^2], but with only the r]/i] — > 77 
modes. In the present framework extended to the 77/77' — > 7r + 7r~7 decay modes, the 
conclusion is confirmed : The recommended value is still clearly preferred ; fit quality 
indicates that it could be slightly smaller (in the direction of the Primakoff measurement), 
but larger values (in direction of the 77 measurements) are clearly disfavored. 



12 Summary And Conclusion 

The conclusions we get are of various kinds. Therefore, we prefer segmenting into 
Subsections. 

Experimental Relevance Of The Box Anomaly 

Concerning the analysis of a possible occurence of the box anomaly phenomenon in 
77/77' decays, the main results reported in the present paper can be summarized as follows : 

• There is a strong evidence in favor of a contact term contribution in the 77/77' 
decays to 7r + 7r _ 7. All aspects (invariant-mass spectra and partial widths) of the 
77/77' — ► 7r + 7r~7 decays can be predicted with fair accuracy using a few pieces of 
information coming from fits to VP7 and Ve + e~ decay modes in isolation and 
from information coming from fit to the pion form factor. 

• The needed contact term is numerically at the precise value predicted for the box 
anomaly contribution by the anomalous HLS Lagrangian. This plays a crucial 
role in yielding, without any fit, the correct dipion invariant-mass spectra and the 
correct partial widths for both the 77 and 77' mesons. 

• If one lets free the parameters defining the p meson lineshape in the 77/77' spectra, 
they stay very close to the values expected from (independent) fits to the pion form 
factor if the predicted contact term is switched on. 



29 



In contrast, if one removes this from the amplitudes, the decription is poor and 
can only be improved by letting the p lineshape becoming inconsistent with what 
is expected from fits to the pion form factor. 

• The fit value obtained for the single free parameter (x, accounting essentially for 
nonet symmetry breaking) indicates indubitably that a global description of all 
VP r y modes and of the four 77/77' decay modes examined here is derived with no 
additional free parameter. 

This leads us to conclude to a strong evidence in favor of the occurence of the box 
anomaly phenomenon in 77/77' — > 7r + 7r~7 decays at precisely the level expected from the 
HLS Model and the WZW Lagrangian. 

Anomaly Equations And Mixing Angles 

On the other hand, we have been led to reexamine the validity of the one-angle 
traditional equations giving the amplitudes for 77/77' — > 77 and — > 7r + 7r~7 at the 
chiral point, when breaking flavor SU(3) and Nonet Symmetries ; resp. the triangle and 
box anomaly equations. 

We have found that the broken HLS Model leads to one-angle (Op) expressions for the 
anomaly equations which match low energy QCD expectations as expressed by (E)ChPT, 
but are in deep contradiction with the equations traditionally used. 

Instead of depending on three unconstrained parameters, this set of (four) equations 
we get depends on only one parameter, closely associated with Nonet Symmetry breaking 
(called a; or A in the body of the text) ; they also depend on z = [fx/ fw} 2 which can 
hardly be considered as a free parameter. They are proved to be easily fulfilled by the 
relevant 77/77' partial widths with fair accuracy. 

Relying on the condition Oq = 0, well accepted by the existing data, the broken HLS 
Model leads to an expression of Op in terms of z and x (or A) which can be approximated 
by a simple formula. Additionally, under the same assumption, an equation leading to 
# 8 ~ 29 P can be derived. 

These equations have been derived from within the framework of the Hidden Local 
Symmetry Model appropriately broken. The phenomenological success of this mecha- 
nism implies that the BKY SU(3) Symmetry breaking scheme, supplemented with Nonet 
Symmetry breaking can be considered as the relevant breaking mechanism. 

This extended BKY breaking scheme forces to a field transformation which admits a 
reliable approximation valid at leading order in the breaking parameters ([z — 1], [x — 1]). 
The refinements permitted by the exact transformation are found beyond the present 
accuracy of the experimental data. 

Perspectives 

At the level of accuracy permitted by the existing data, the HLS Model (including its 
anomalous sector) together with the extended BKY symmetry breaking scheme, covers 
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successfully all aspects of the experimental data examined so far, certainly up to the 
mass. 

It would be interesting to have improved data in order to check up to which accuracy 
the HLS framework is predictive. For this purpose, more and better data on the 77/77' 
sector would be welcome. 

These could come from tau-charm factories (CLEO-C and upgraded BESS) which, 
running at the J/tp(lS), produce very large samples of 77/77' mesons under especially clean 
physics conditions. For instance, in the run at the J ftp (IS) foreseen by CLEO-C in 2005 
10 9 events will be collected. This will provide 860 000 77 produced opposite in azimuth 
to a single monoenergetic photon and 2 000 000 opposite to w/(f>. The corresponding 7/ 
decay modes will provide samples of about 4 300 000 77' produced opposite to a single 
photon and 500 000 opposite to u/<j>. This should allow an exhaustive study of the 77/77' 
system and a much better understanding of low energy QCD. 
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Appendix A 



Al : The Exact Field Transformation 



As stated in Section the field transformation given by Eq. (JSJ) is an approximation 
of the full transformation which has been derived in [32J. 

In order to bring the kinetic energy part of the U(3)/SU(3) broken HLS Lagrangian 
into canonical form (see Eq. (JJJ)), it is appropriate to perform the renormalization in two 
steps. One first diagonalizes the standard Chls piece using the field transformation Eq. 
(JHJ). This makes the Lagrangian canonical for the tt/K sector, and one yields intermediate 
fields for the isoscalar sector (double prime fields). In terms of bare fields, we have : 



zr 



cos (3 
- sin (3 



sin (3 



cos (3 y= sin (3 

v2 



(A.l 



where, one has defined : 



(2^ + l) 2 + 2(z-l) 
3z 



0.90 



tan [3 = V2 



2z 



0.20 



This transformation brings the kinetic term in the following form 

2 T = [dr)l\ 2 + [<9r^] 2 + Ar [sin/? dq'i + cos/? d^} 2 . 



(A . 2 



(A . 3 



The transformation to fully renormalized fields (primed fields) is performed with : 

(A. 4) 

where v carries the real information about Nonet Symmetry breaking (see Eq. ( |7| l) : 



V's' 




1 + v sin 2 (3 


v sin p cos (3 






Vo . 




v sin (3 cos (3 


1 + v cos 2 (3 







VI + Ar 2 - 1 ~ 0.10 



(A . 5 



That the transformation combining Eqs. (|A . 1 j) and (|A . 4 \j results, at leading order 
in the breaking parameters [z — 1] and [x — 1], into a transformation as simple as Eq. 
(JBJ), is a little bit unexpected. As noted in the main text, there are several combinations 
involving A which are equivalent to x at leading order ; they are all of the form exhibited 
by Eq. (0) which is typically a good representation of x in terms of A. 

This remainder makes clear why x is influenced by the SU(3) symmetry breaking. A 
typical expression for x is : 

x = . 1 (A . 6 ) 
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A2 : The Anomalous Amplitudes At The Chiral Limit 



The expressions for the anomalous amplitudes at the chiral limit, when using the 
exact transformation, are easy to get. They amount to the following changes for the 
triangle anomaly expressions in Eqs. : 



hz-2 
3z 

v2— — x 

6Z 



1 



1+V 
1 

l + « 



2 v cos 3 
+ - 



V2 



3z rz 
5z + 1 v sin 3 



3 2 



rz 



(octet) 



(singlet) 



(A. 7) 



The octet and singlet combinations for the box anomalies can easily be identified by 
the occurence of the x factor in Eqs. (|T9"j) and (|21|). The changes to be performed there 
are : 



1 



x 



1 



1 



1 



1 + 



1 - 



v cos (5 



rz 



v sin 3 



(octet) 



(singlet) 



(A . 8 



rz^/2 

It is worth remarking that the exact field transformation changes the (pjrj) and (pjrf) 
coupling constants in such a way that the ex's -modified as just stated- still factor 
out from their expression. Therefore, Eqs. ()18|) and (J23j) keep their structure and the 
decay invariant-mass spectra for the rj/i]' are the same as for the approximate field 
transformation. 



A3 : Decay Constants And Mixing Angles 

One can express easily the EChPT coupling constants (/o and /§) and mixing angles 
(6*o and 9$) in terms of the parameters mixing A and 3 defined in the previous Subsections. 

The following matrix elements of axials currents can be defined in the broken HLS 
Lagrangian [32] : 

(0J> 8 (g)) = z/ 8 g M , (0|J>°(?)> = i/o^ 

(0|J>o(g)> = *b 8qfl , (0| J> 8 (g)> = ib q„ (A . 9 ) 

One can easily write down the currents and their matrix elements (see |32]> Section 6) 
in the case when the field transformation is not approximated by Eq. (jSJ). Using the 
notations defined in [32], one finds first : 

rz 

[1 + v cos 23 ] cos (3 

1 + v 

(A . 10 ) 

TZ 

■ [1 — v cos 2f3 ] sin f3 



f* 
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Defining the following parameter combinations 



hi = cos (3 — 



one also yields 



sin (3 
~V2 



0.90 



A cos (3 — —7=- sin (3 ~ 0.03 
v2 



rz 
l+v 



hi (1 + A) + v(2 cos (3 + h 2 ) sin 2 (3 



fo 

u 

■j- = — sin/3- [ (1 + t>cos2/3) — /iit>cos/3 ] 



(A . 11 



(A . 12 



A4 : The Condition 9a = 



Phenomenology [oz] as well as explicit EChPT computations [58] indicate that the 
mixing angle #0 is very close to zero. Table El clearly illustrates that present data are 
statistically insensitive to letting 80 departing from zero. Under such conditions, several 
interesting relations show up. 

The definition of the angles #0 and can be expressed in terms of the parameters in 
Eq. ()A . 9 jl [32]. Using Eq. (|A .12 |) . one can derive : 



tan #8 = tan (6p + ip$) 



tan 8q = — tan {dp — tpo) 



(A . 13 



where tanyj 8 = b$/ f$ and tan<^ = &0//0 can be explicitly computed. The condition 
#o = strictly implies that 9p = ifo which gives : 



tan#/ 



1 + A 



tan f3 



^tan/3 



v tan (3 
1 + fJ- + 



(A . 14 



From Eq. ()A . 6 |) . the first term can be interpreted as x 2 and the product of the 
first two factors is just Eq. (jlO)) for tan#p modified with x 2 . With the values for v 
and tan/3 we have mentioned, the leading correction amounts to only 1.5 10~ 2 . If one 
keeps a l/x/l + A (corresponding to having x in Eq. ()10[0 in front of this expression, the 
correction term gets a additional contribution —A/2 ~ 5. 10 -2 which becomes dominant. 
Therefore : 



tan 9 1 



V2- 



1 



2 + z 



X' 



(A . 15 



could indeed be prefered to Eq. ([TU]). 

The second information which follows from 6 = is an approximate relation between 
6g and the wave-function mixing angle 8p : 



tan 8s = 2tan#p 
where the leading correction is ~ 7. 10~ 2 . 



tan (3 
2V2 



+ ■ 



(A . 16 
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Layter [TT] 





ei 


e 2 

(GeV- 2 ) 


X 2 /dof 
(Prob.) 


Tj 

PW. (eV) 


PW. (keV) 


CT + No Fit 


0.222 ± 0.011 


-1.203 ± 0.017 


37.88/35 
(34%) 


56.3 


48.9 


CT + Fit 


u.ooa_ 056 


i.oyo_ . 160 


35.16/33 
(36.6%) 


53.3 


46.2 


No CT + No Fit 


0.222 ±0.011 


-1.203 ±0.017 


140.65/35 
(0%) 


100.9 


57.5 


No CT + Fit 


u.yoo_ 093 


9 qc;c;+o.180 
z,.ooo_ Q 200 


60.88/33 
(0.2 %) 


75.3 


40.0 


Gormley |10j 




ei 


e 2 

(GeV- 2 ) 


X 2 /dof 
(Prob.) 


V 

PW. (eV) 


i 

PW. (keV) 


CT + No Fit 


0.222 ±0.011 


-1.203 ±0.017 


25.58/34 
(85%) 


56.3 


48.9 


CT + Fit 


0.269 ±0.080 


i 97c -+0.135 
l.Z / O_o.l55 


25.01/32 
(80.6 %) 


54.6 


47.7 


No CT + No Fit 


0.222 ±0.011 


-1.203 ±0.017 


54.13/34 
(1.6%) 


76.9 


53.4 


No CT + Fit 


0.529 ±0.090 


-1 700+ - 154 

L - 1 uu -0.195 


36.38/32 
(27.2%) 


65.7 


44.5 



Table 4: Simultaneous fits of the 77/77' distributions from |2l^lTT] on the one hand, and 
from [2*3*1 ITU] on the other hand. CT stand for the "contact terms" generated by the box 
part of the WZW Lagrangian (see Eq.(jSJ)). The values for ei and e 2 quoted in "No Fit" 
entries are taken from Table d and not varied from their central values. The " PW." 
entries are the central values for the 77 and 77' partial widths in the appropriate units ; 
the recommended values for these |Hj are given in Table El 
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PDG 2002 Fit Result Significance (n a) 

x 0.911 ±0.015 

X 2 /dof 2.66/3 

Probability 44.6 % 





•> 77) (keV) 


0.46 ±0.04 


0.46 ±0.01 


0.00 a 


T(rf- 


-> 77) (keV) 


4.29 ±0.15 


4.34 ±0.14 


0.24 a 


T( V - 


■> 7r + 7r _ 7) (eV) 


55 ±5 


56.64 ± 1.71 


0.31 a 


T(rf- 


-> 7r + 7r _ 7) (keV) 


60 ±5 


49.75 ±3.88 


1.62 (7 



Table 5: Simultaneous fit of the four HLS anomaly equations (Eqs. ([11))) and ([21)1 with 
only x free (Approximate Field Transformation). First data column gives the recom- 
mended values 0. 







PDG 2002 


Fit Result 
#0 free 


Fit Result 
#o = 


A 






0.23 ±0.06 


0.21 ±0.04 


o P 






-10.85° ± 1.27° 


-11.48° ±0.02° 


X 2 /dof 




2.93/2 


3.20/3 


Probability 




23.1 % 


36.5% 








-1.01° ± 1.27° 











-19.37° ± 1.29° 


-18.16° ±0.24° 


/o 






1.37 ±0.03 


1.36 ±0.03 


fs 






1.34 ±0.01 


1.34 ±0.02 




77) (keV) 


0.46 ± 0.04 


0.45 ±0.03 


0.44 ±0.01 




77) (keV) 


4.29 ±0.15 


4.37 ±0.23 


4.20 ±0.17 


rfa- 


7r + 7r~7) (eV) 


55 ±5 


55.38 ±2.78 


55.98 ± 1.76 


r(i/- 


7r + 7r _ 7) (keV) 


60 ±5 


49.40 ±4.85 


46.93 ±4.00 



Table 6: Simultaneous fit of the four HLS anomaly equations modified by using the exact 
field transformation. The second data column reports on letting free A and dp, while in 
the third data column only A is allowed to vary. 
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Arbitrary Units 



T}' No Contact Term 
rf Contact Term 




0.8 
0.6 
0.4 
0.2 




Arbitrary Units 



c) 




7] No Contact Term 
■f] Contact Term rescaled 



0.2 



0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 

m„„ GeV/ c 2 



Figure 1: 

Predicted shapes for rf (top) and r\ (mid) distributions as functions of the dipion 
invariant mass. Full line histograms correspond to having the contact term in the 
amplitude, dotted line histograms correspond to removing the contact term from 
the amplitude. All other numerical parameters are at the same values (see Table [Q). 
In the bottom figure, we plot the prediction when accounting for the contact term 
(rescaled) superimposed with the prediction derived by removing this contribution. 
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Figure 2: 

Invariant dipion mass Distributions for rj' decay. Experimental data 
sets with the predicted distributions without the contact term (dashed curve) 
and with this contribution activated (full curve). The numbers given are 
X 2 (contact term) /x 2 (no contact term) for the lineshapes only. 
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Number of events 



Number of events 




Figure 3: 

Invariant dipion mass Distributions for rj' decay. Experimental data 
sets with the predicted distributions without the contact term (dashed curve) 
and with this contribution activated (full curve). The numbers given are 
X 2 (contact term) /x 2 (no contact term) for the lineshapes only. 
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Figure 4: 

Invariant dipion mass Distributions for rj' decay and the single rj decay 
(as a function of the photon momentum in the 77 rest frame). Experimen- 
tal data sets with the predicted distributions without the contact term (dashed 
curve) and with this contribution activated (full curve). The numbers given are 
X 2 (contact term) /x 2 (no contact term) for the lineshapes only. 



45 



Number of events 



2250 
2000 
1 750 
1500 
1250 
1000 
750 
500 
250 




700 
600 
500 
400 
300 
200 
1 00 




• r, (LAYTER) 

p- 



7] Contact Term 
7] No Contact Term 



0.06 



0.08 



0.12 



0.14 



Arbitray Units 



0.18 0. 

K v GeV/c 2 



~- • r, (GORMLEY) 
(5.5/12.3) 



rj Contact Term 
v No Contact Term 



0.06 



0.08 



0.1 



0.12 



0.1 4 



0.' 6 



0. ' 8 

K, GeV/c 2 



Figure 5: 

Photon momentum distribution in 77 decay Experimental data are from 
Layer et al. jTT] (top), and Gormley et al. [TU] (bottom) ; Experimental data 
sets with the predicted distributions without the contact term (dashed curve) 
and with this contribution activated (full curve). The numbers given are 
X 2 (contact term) /x 2 (no contact term) for the lineshapes only. 
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Figure 6: 

Fit of the 7/ invariant mass spectrum, top by including the contact term, bottom 
by removing this term. Compare the peak locations. 
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